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Abstract. Quantum chaos is a subject whose major goal is to identify and to investigate different 
quantum signatures of classical chaos. Here we study entanglement production in coupled chaotic 
systems as a possible quantum indicator of classical chaos. We use coupled kicked tops as a model for 
our extensive numerical studies. We find that, in general, presence of chaos in the system produces 
more entanglement. However, coupling strength between two subsystems is also very important 
parameter for the entanglement production. Here we show how chaos can lead to large entanglement 
which is universal and describable by random matrix theory (RMT). We also explain entanglement 
production in coupled strongly chaotic systems by deriving a formula based on RMT. This formula 
is valid for arbitrary coupling strengths, as well as for sufficiently long time. Here we investigate 
also the effect of chaos on the entanglement production for the mixed initial state. We find that 
many properties of the mixed state entanglement production are qualitatively similar to the pure 
state entanglement production. We however still lack an analytical understanding of the mixed state 
entanglement production in chaotic systems. 
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1. Introduction 

Entanglement is a unique quantum phenomenon which can be observed in a system con- 
sists of at least two subsystems. In case of an entangled system even if we know the exact 
state of the system, it is not possible to assign any pure state to the subsystems and that 
leads to the weU-known unique quantum correlations which exists even in spatially well 
separated pairs of subsystems. This phenomenon was first discussed by Schrodinger to 
point out the nonclassicality implied by the quantum mechanical laws [1]. This remark- 
able feature of quantum mechanics has recently been identified as a resource in many areas 
of quantum information theory including quantum teleportation [2], superdense coding [3] 
and quantum key distribution [4]. Moreover, entanglement is also a key ingredient of all 
the proposed quantum algorithms which outperform their classical counterparts [5,6]. 

Quantum mechanical study of classically chaotic systems is the subject matter of 'quan- 
tum chaos' [7,8]. A major challenge of quantum chaos is to identify quantum signatures of 
classical chaos. Various signatures have been identified, such as the spectral properties of 
the generating Hamiltonian [9], phase space scaiiing [10], hypersensitivity to perturbation 
[11], and fidelity decay [12], which indicate presence of chaos in underlying classical sys- 
tem. Recent studies have shown that entanglement in chaotic systems can also be a good 
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indicator of the regular to chaotic transition in its classical counterpart [13-24]. A study 
of the connections between chaos and entanglement is interesting because the two phe- 
nomena are prima facie uniquely classical and quantum, respectively. This is definitely an 
important reason to study entanglement in chaotic systems. Moreover, presence of chaos 
has also been identified in some realistic model of quantum computers [25]. 

In this paper, we have investigated entanglement production in coupled chaotic systems. 
We have used coupled kicked tops as a model for our whole study. We have considered 
the entanglement production for both chaotic and regular cases. Moreover, we have also 
considered the effect of different coupling strengths on entanglement production. Most of 
the earlier studies have considered the effect of chaos on entanglement production for the 
case of initially pure state of the overall system. A basic assumption of these studies is 
that the initial state of the overall system is completely known. However, in many of the 
realistic scenarios, we do not have a complete knowledge of the state of a quantum system. 
For instance, when a quantum system interacts with its surroundings, it is not possible 
to know the exact state of the system. We may only express the state of the system as 
statistical mixture of different pure states, and that is a mixed state. In this paper we have 
also studied mixed state entanglement production in chaotic systems. 

This paper is organized as follows. In the next section we discuss about classical and 
quantum properties of two coupled kicked tops, our primary model. Then we have defined 
the measures of both pure and mixed state entanglement. Finally, we have concluded this 
section with a discussion on the initial states (both pure and mixed) used here. In Sec. 3., 
we present the numerical results on the entanglement production in coupled kicked tops 
for different single top dynamics and also for different coupling strengths. Here we have 
studied entanglement production for both pure and mixed initial state. In Sec.4., we derive 
the statistical universal bound on entanglement using random matrix theory (RMT). We 
also derive an approximate formula, based on RMT, to explain the entanglement production 
in coupled strongly chaotic systems. Finally, we sunraiarize in Sec.5. 

2. Preliminaries 

2. 1 Coupled kicked tops 
1 Quantum top 

The single kicked top is characterized by an angular momentum vector J = {J^, Jy, Jz), 
where these components obey the usual conamutation rules. The Hamiltonian of the single 
top is given by [26] 



The first term describes free precession of the top around y axis with angular frequency 
tt/2, and the second term is due to periodic (5-function kicks. The second term is torsion 
about z axis by an angle proportional to J^, and the proportionality factor is a dimen- 
sionless constant k/2j. Now the Hamiltonian of the coupled kicked tops can be written, 
following Ref. [14], as 




n=—OG 



(1) 



2 



n{t) = Hi{t) + H2{t) + Huit), (2) 

where 

Hi2 (t) = -. J,, J., V <5(i - n) , (3b) 

wherei = 1,2. Here iJi(t)'s represent the Hamiltonians of the individual tops, and -ffi2(i) 
is the couphng between the tops via spin-spin interaction with a couphng strength e/j. 
Corresponding time evolution operator, defined in between two consecutive kicks, is given 
by 

Ut = Ul^iUi U2) = Ul^iiuM) ® iPM)], (4) 
where the different terms are given by, 

[// = exp (^-i^Jy?) , Uf = exp > ^12 = exp [-i-JziJz^ (5) 

and as usual i = 1,2. 



2 Classical top 

The classical map corresponding to the coupled kicked tops can be obtained from the quan- 
tum description with the Heisenberg picture in which the angular momentum operators 
evolve as 

J„+i = u\.3JJt. (6) 

ExpUcit form of this angular momentum evolution equation for each component of the 
angular momentum is presented in Ref. [17]. We now proceed by rescaling the angular 
momentum operator as (Xj, Y^, Zi) = {Jxi, Jyi, Jzi)/ h for i ~ 1, 2. The commutation 
relations satisfied by the components of this rescaled angular momentum vector as follow 
: [Xi, Yi] = iZi/j, \Yi, Zi] = iXi/j and [Zj, Xi] = iYi/j. Therefore, in j ^ 00 hmit, 
components of this rescaled angular momentum vector will commute and become classical 
c-number variables. In this large-j hmit, we obtain the classical map corresponding to 
coupled kicked top as [17]: 



X[ = Zi cos A12 + Yi sin A12, (7a) 

Yl = -Zi sin A12 + Fi cos A12 , (7b) 

Z[ = -Xi, (7c) 

X'2 = Z2 cos A21 -I- Y2 sin A21 , (7d) 

F2 =~Z2 sin A21 + Y2 cos A21 , (7e) 

Z'2 = -X2, (7f) 
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Figure 1. Phase space pictures of the single top, corresponding to different parameter 
values, are presented, (a) k = 1.0. Phase space is mostly covered by the regular region, 
(b) k = 2.0. The phase space is still very much regular, but now a thin stochastic layer 
can be observed at the separatrix. (c) k = 3.0. The phase space is truly mixed type. 
Few regular elliptic islands are visible inside the chaotic region, (d) k = 6.0. The 
phase space is almost covered by the chaotic region with few tiny elliptic islands. The 
soUd circle (•) is the point at which we will construct the initial wave packet during our 
study of the pure state entanglement production 



where 

Ai2 = kXi + eX2 and A21 = kX2 + eXi. (8) 

In the hmit e ^ 0, the classical map for the coupled kicked tops decouple into the classical 
map for two single tops. The classical map for one such uncoupled top can be written as 

X' = Z cos kX + Y sin kX (9a) 
Y' = -ZsmkX + Y coskX (9b) 
Z' = -X. (9c) 

It is clear from the above expression that the variables {X, Y, Z) lie on the sphere of ra- 
dius unity, i.e. X^ + Y"^ + Z'^ = 1. Consequently, we can parameterize the dynamical 
variables in terms of the polar angle 9 and the azimuthal angle as X = sin cos </>, F = 
sin 6 sin (f), and Z = cos 9. In Fig. 1, we have presented the phase space diagrams of the 
single top for different values of the parameter k. For k = 1.0, as shown in Fig. 1(a), the 
phase space is mostly covered by regular orbits, without any visible stochastic region. Our 
initial wave packet, marked by a solid circle at the coordinate (0.89, 0.63), is on the reg- 
ular elUptic orbits. As we further increase the parameter, regular region becomes smaller. 



Fig. 1(b) is showing the phase space for k = 2.0. Still the phase space is mostly covered 
by the regular region, but now we can observe a thin stochastic layer at the separatrix. In 
this case, the initial wave packet is on the separatrix. For the change in the parameter value 
from = 2.0 to = 3.0, there is significant change in the phase space. At A; = 3.0, shown 
in Fig. 1(c), the phase space is of a truly mixed type. The size of the chaotic region is now 
very large with few regular islands. At this parameter value, the initial wave packet is in- 
side the chaotic region. Fig. 1(d) is showing the phase space for k = 6.0. Now the phase 
space is mostly covered by the chaotic region, with very tiny regular islands. Naturally, 
our initial wave packet is in the chaotic region. 

2.2 Measures of entanglement 
1 Pure state 

Entanglement measure for a system consisting of two subsystems (bipartite) is well defined 
if overall state of the system is in a pure state. In this case subsystem von Neumann 
entropy, i.e. von Neumann entropy of the reduced density matrices (RDMs), is a measure 
of entanglement. If there is no entanglement among the two subsystems, then the RDMs 
will correspond to density matrices of pure states and hence the subsystem von Neumann 
entropy will vanish. Otherwise, in case of entanglement, a non-zero value of the subsystem 
von Neumann entropy will be a measure of entanglement among the two subsystems. 

Let us assume that the state space of a bipartite quantum system isTi = Hi Ti.2, 
where dim Hi = < dim 7^2 = M, and dimH = d ^ NM. If p = Ei-Pil'/'* 
is an ensemble representation of an arbitrary state in H, the entanglement of formation 
is found by minimizing J^iPi-^ il'Pi)) o^^'" possible ensemble realizations. Here E 
is the von Neumann entropy of the RDM of the state \(pi) belonging to the ensemble, 
i.e. its entanglement. For pure states \tp) there is only one unique term in the ensemble 
representation and the entanglement of formation is simply the von Neumann entropy of 
the RDM. 

The two RDMs of any bipartite pure state \ip) are pi = TT2{\ip){tp\) and p2 = 
Tri{\tp){^p\). The Schmidt decomposition of \tp) is the optimal representation in terms 
of a product basis and is given by 

N 

|V') = ^ VA.|</-(^))|0f), (10) 
j=l 

where < Aj < 1 are the (non-zero) eigenvalues of either RDMs and the vectors are the 
corresponding eigenvectors. The von Neumann entropy Sy is the entanglement E{\ip)) is 
given by 

N 

Sv = -Tri{pilnpi) = -J2Xiln{Xi) ■ 1 = 1,2. (11) 

The von Neumann entropy can only be calculated in the eigenbasis of the RDMs due to 
the presence of logarithmic function in its definition. Therefore it is not easy to calculate 
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this measure unless one has some information of the eigenvalues of the RDMs. Conse- 
quently linearized version of the von Neumann entropy, called linear entropy, has also 
become a popular measure of entanglement. This measure of entanglement is defined as 

N 

Sr = 1- TxipI = 1 - ^ ; Z = 1,2. (12) 

i=l 

The linear entropy can be calculated without knowing the eigenvalues of the RDMs, be- 
cause Trjp^ is equal to the summation of absolute square of all the elements of RDMs. 
However, strictly speaking, the linear entropy is not a true measure of entanglement, rather 
it is a measure of mixedness of the subsystems which increases with entanglement among 
the two subsystems. Therefore, the linear entropy can be considered as an approximate 
measure of entanglement. 

2 Mixed state 

A major issue related to the study of the mixed state entanglement is lack of unique mea- 
sure of entanglement. Probably this issue has discouraged any work related to the mixed 
state entanglement production in chaotic systems. Recently Vidal and Werner [27] have 
proposed a computable measure of entanglement called Log-negativity following Peres' 
criterion of separability [28]. We use this measure to characterize mixed state entangle- 
ment production in chaotic systems. Basic idea of this measure is very simple and straight- 
forward to state. 

A most general form of a separable bipartite mixed state is given by 

p=E«'^i'^®'^^^ (13) 

i 

where the positive weight factors Pi satisfy ^ ■ pi = 1, p\^^ and p\^^ are density matrices 
for the two subsystems. We can construct a matrix p^^ from p by taking transpose only 
over second subspace, i.e. 




i 



This partial transpose operation is definitely not a unitary operation, but is still Hermi- 
tian. The transposed matrices (^P^^^ ^re positive matrices, and hence they are legitimate 

density matrices. Consequently, if p is separable, p'^^ is a positive matrix. This is also true 
for . In general, this is a necessary condition of separability. 

Log-negativity measures the degree to which p^^ (or p^^ ) fails to be positive. If p is 
an entangled state, then may have some negative eigenvalues. The Log-negativity is 
logarithm of the sum of absolute value of the negative eigenvalues of p^^ which vanishes 
for unentangled state. It can be shown by simple algebraic manipulation that the sum 
of absolute value of all the negative eigenvalues of p^^ is linearly related to the sum of 
absolute value of all the eigenvalues of p^'^ . Therefore, the Log-negativity measure [p) 
can be defined as 
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ENip) = \n(j2\K\^ (15) 
where d is the dimension of p. 

2.3 Initial state 
1 pure state 

We use generahzed SU(2) coherent state or the directed angular momentum state [7,26] as 
our initial state for the individual tops and this state is given in standard angular momentum 
basis \ j, m) as 




{3,m%Ao) = {l + \lV)-'l'-'^\ [ ,7' K (16) 



where7 = exp(i^o) tan(^o/2). For the coupled kicked top, we take the initial state as the 
tensor product of the directed angular momentum state corresponding to individual top, 
i.e., 

|V'(0)) = |^W,<^^^^)|^f ,4^)), (17) 

where (^q*-* , (/)q*'' ) = (0.89,0.63) for i = 1.2. This initial state is evolved under Ut as 
\'ip{n)) = U^\'tp(Q)) for different values of the parameter k and for different couphng 
strength e, and the results are displayed in Fig. 2. 



2 mixed state 



In this case we have considered a very simple unentangled mixed state, where the initial 
state corresponding to first top is mixed and the same corresponding to the second top is 
pure. Mathematically we express this state as p(0) = Pi(0) ® |V'2(0)) (■i/'2(0)|, where 
PiiQ) is the initial mixed state of the first subsystem and \ip2{Q)) is the initial pure state of 
the second subsystem. We take |V'2(0)) as a generalized SU(2) coherent state as presented 
above in Eq.(16). The mixed state /5i(0) is a combination of two such coherent states 
placed at two different points on the phase space, i.e., 

pi(o) =p|^?«,0«) + (i-p) \e'iS.<t>'-^) (18) 

Here we choose (6'io\tfo) = (0.89, 0.63) and {e\^^Am) = (2.25, -0.63) in such a way 
that the dynamical properties of these points are similar for any value of k. For the second 
top, we choose |V'2(0)) = |6'^^' = 0.89,4^' = 0.63). We only consider p = 1/2 case, 
this means the contribution of each coherent state is same on the formation of pi{Q)- The 
initial mixed state p(0) is evolved under Ut as p(n) = U^p{0)U^". We study the time- 
evolution of the Log-negativity measure for different k and e, and the results are displayed 
in Fig. 3. 
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3. Numerical results 

3.1 Pure state entanglement production 

In Fig.2, we have presented our results for the entanglement production in coupled kicked 
tops for the spin j = 80. As we go from top to bottom window, coupling strength is 
decreasing by a factor of ten. Top window corresponds to e = 10^^, middle window is 
showing the results for e = 10~^, and the bottom one corresponds to the case e = 10"''. 
For each coupUng strengths, we have studied entanglement production for four different 
single top parameter values, whose corresponding classical phase space picture has already 
been shown in Fig.l. 

1 Coupling e — 10~^ 

The entanglement production for this strong coupling strength has been presented in 
Fig.2(a). It shows that there exists a saturation of Sv for the regular cases (fc = 1.0 
and k = 2.0), which are much less than the saturation value corresponding to strongly 
chaotic cases such as when k ~ 6.0. The saturation value of Sy for k = 6.0 is a statistical 
bound Sv = In(iV) ~ \ ^ 4.57 (where N = 2j + 1 = 161), which can be estimated 
analytically from RMT [16], and we will discuss about this in the next section. However 
for fc = 3.0, corresponding to a mixed classical phase space, the saturation value of Sy 
is less than the above mentioned statistical bounds, which indicates the influence of the 
regular regions. These distinct behaviors of the entanglement saturation can be understood 
from the underlying classical dynamics. For k = 1.0, the initial unentangled state is the 
product of the coherent wave packet placed on some elUptic orbits of each top. Therefore, 
the evolution of this unentangled state under the coupled top unitary operators is restricted 
by those elliptic orbits. Finally, the wave packet spreads all over those elliptic orbits and 
the entanglement production reaches its saturation value. At fc = 2.0, the center of the 
initial coherent state is inside the separatrix. Therefore, in its time evolution, the spreading 
of the wave packet is restricted to be inside the separatrix region. Finally it spread over 
the whole separatrix region, and the entanglement production arrives at its saturation. At 
k = 3.0 and k = 6.0, the initial wave packets are inside the chaotic region. However, 
due to the smaller size of the chaotic region corresponding to the case of fc = 3.0 than the 
case corresponding to fc = 6.0, the wave packet can spread over less of the phase space for 
fc = 3.0 than fc = 6.0. Consequently, the saturation value of the entanglement production 
is less for fc = 3.0 than fc = 6.0. 



2 Coupling e = 10 ^ 

Let us now discuss the case of coupling strength e = lO"'^, whose results are presented 
in Fig.2(b). For the non-chaotic cases (fc = 1.0 and fc = 2.0), the saturation value of 
the entanglement production is less than the entanglement saturation value observed in 
the stronger coupling case (e = 10"^). This is because, for weaker coupling case, the 
interaction between two subsystems is less and the individual subsystems behave more like 
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Figure 2. Time evolution of ttie von Neumann entropy in coupled kicked tops is pre- 
sented for different coupling strengths and for different underlying classical dynamics, 
(a) e = 10"^ (b) e = 10"^ (c) e = 10"*. SoM line represents k = 1.0, dotted 
line corresponds to fe = 2.0, dashed line is for fe = 3.0 and dash-dot line represents 
k = 6.0. 

isolated quantum systems. Similarly, for the strong chaos case {k = 6.0), the entanglement 
production is well short of the known statistical bound ln(A'^) — i. 

3 Coupling €= 10~^ 

The entanglement production for this very weak coupling strength has been presented in 
Fig.2(c). The entanglement production for the weakly coupled strongly chaotic system has 
recently been explained by perturbation theory [18]. However, the formula presented in 
that work is only valid for short time. In the next section we have presented an approxi- 
mate formula for the entanglement production in coupled strongly chaotic systems which 
is valid for sufficiently long time and for any arbitrary coupling strengths. This formula ex- 
plains the entanglement production for the strongly chaotic case {k = 6.0). Here we have 
observed an interesting phenomenon that the entanglement production is much larger for 
the non-chaotic cases than the chaotic cases. Rather, we can say that, for weakly coupled 
cases, the presence of chaos in the systems actually suppresses entanglement production. 



3.2 Mixed state entanglement production 

In Fig. 3, we have presented the Log-negativity measure Em{p) of the mixed state entan- 
glement production for different individual top dynamics [k = 1.0, 2.0, 3.0, and 6.0) and 
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Figure 3. Evolution of the Log-negativity measure (evolving under the coupled tops 
time evolution operator Ut)- Solid lines and dotted lines are representing the results 
corresponding to the non-chaotic cases (fc = 1.0 and k = 2.0, respectively). Dashed 
lines are representing the mixed case (k = 3.0) and dash-dot lines are showing the 
results for the strongly chaotic case (fc = 6.0). (a) Representing the results for the 
stronger coupling strength (e = 1.0). (b) This window is showing the results for 
e = 0.1. (c) This window is for e = 0.01 case, (d) This window is showing the 
results for the weak coupling case (e = 0.001). 

for different coupling strengths. 



1 Coupling e = 1.0 

Let us start the discussion with the case of strong coupling e = 1.0, whose results are pre- 
sented in Fig. 3(a). This couphng strength is so strong that, irrespective of the individual top 
dynamics, the overall coupled system is chaotic. Therefore, the location of the initial state 
and the dynamics of the individual tops are irrelevant for the saturation of Em{p)- Con- 
sequently, we have observed almost same saturation value of Em{p) for all the different 
individual top dynamics. 

2 Coupling e = 0.1 

The time evolution of Em{p) corresponding to e = 0.1 is presented in Fig. 3(b). For 
this coupling strength, we have observed that the saturation value of Em (p) for the non- 
chaotic cases {k = 1.0 and A; = 2.0) are less than the saturation value corresponding to 
other two cases. These lower saturation values of En {p) for the non-chaotic cases indicate 
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the influence of the regular orbits on the mixed state entanglement production. We have 
also noticed for the non-chaotic cases that the saturation value of E^ip) is less than the 
saturation value observed in the stronger coupling case (e = 1.0). However, the saturation 
value of En (p) corresponding to other two cases, k = 3.0 and k = 6.0, are almost equal 
to the previous case (e = 1.0). 

3 Coupling e= 10~^ 

The mixed state entanglement production for the coupling strength e = 10~^ has been 
presented in Fig. 3(c). Here again the saturation value of En{p) corresponding to the non- 
chaotic cases are less than the other two cases. Moreover, due to the weaker coupling, 
the saturation value of En{p) for the non-chaotic cases are less than the saturation value 
observed in the previous two cases of stronger coupling strengths (e = 1.0 and e = 0.1). 
Here we have first time observed an interesting phenomenon that the magnitude of Em{p) 
corresponding to the non-chaotic cases are larger than the other two cases at least within 
a time interval 50 < n < 100. However, the saturation value corresponding to the case 
of mixed phase space (fc = 3.0) is still less than the saturation value corresponding to the 
chaotic case [k = 6.0). 

4 Coupling e = 10~^ 

Finally, for the weak coupling strength e = 10~^, we have observed completely different 
behaviors of the evolution of En{p) and these results are presented in Fig. 3(d). For this 
coupling strength, E^ip) corresponding to the non-chaotic cases are always higher than 
the chaotic cases within our time of observation. These results imply that, in case of weak 
coupling, the presence of chaos in the system actually suppresses the entanglement pro- 
duction. This suppression of the entanglement production by chaos for the weak coupling 
case was also observed for the pure state entanglement production, which we have already 
discussed in the previous section. 

4. Some analytical results 

4.1 Random matrix estimation of the statistical bound on entanglement 

In our numerical study of pure state entanglement production, we have observed a statis- 
tical bound on entanglement for the strongly coupled strongly chaotic top. An identical 
property has also been observed for the stationary states of coupled standard map [15] 
which indicates its universality. The parameter values considered in our numerical study, 
the nearest neighbor spacing distribution of the eigenangles of Ut is Wigner distributed, 
which is typical of of any quantized chaotic systems [7,8]. Therefore, it is quite reasonable 
to expect that the statistical bound on entanglement can be estimated by random matrix 
modeling. The two RDMs, corresponding to two subsystems, have the structure A'^ A and 
AA"^ , where ^ is a rectangular N x M matrix containing the vector components of the 
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Figure 4. Distribution of tlie eigenvalues of the RDMs of coupled kicked tops, aver- 
aged over all the eigenstates (A'^ = 2ji + 1 = 33). Solid curves correspond to the 
theoretical distribution function Eq.{19). 



bipartite system. We have pointed out the distribution of the eigenvalues of these RDMs as 
[16] 



/(A) = 



max — A) (A — X mm ) 

27r A 
1 / 1 2 



TV 



« = -(l + -±— 1, (19) 



where A e [Amin, Amax], Q = M/N, and Nf{X) dX is the number of eigenvalues within 
A to A + dX. This has been derived under the assumption that both M and N are large. 
Note that this predicts a range of eigenvalues for the RDMs that are of the order of 
For Q 1, the eigenvalues of the RDMs are bounded away from the origin, while for 
Q = 1 there is a divergence at the origin. All of these predictions are seen to be borne out 
in numerical work with coupled tops. 

Fig.4 shows how well the above formula fits the eigenvalue distribution of reduced den- 
sity matrices corresponding to the eigenstates of the coupled tops. Time evolving states 
also have the same distribution. This figure also shows that the probabifity of getting 
an eigenvalue outside the range [Ajjjjji, Amax] is indeed very small. The sum in Sy [see 
Eq.(l 1)] can be replaced by an integral over the density /(A): 



N / /(A)AlnAdA = \n{-fN) = In(iV) + ln(7). (20) 



•^rmn 



The integral in 7 can be evaluated to a generalized hypergeometric function and the final 
result is 
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When the Hilbert space dimension of the subsystems are equal, that is Q = 1, the above 
expression gives 7 = exp(— |), and therefore the corresponding Sy is ln(7V) — ^. This 
is the statistical bound on entanglement found in our numerical study. In another extreme 
case, when M ^ N, that is Q ^ 1> then 7 ~ 1 and the corresponding Sy is equal 
to its maximum possible value [ln(A/^)]. Therefore, the analytical formulation based on 
RMT is able to explain the saturation behavior of quantum entanglement production very 
accurately. 



4.2 Entanglement production in coupled strongly chaotic system 



We have already mentioned that, due to the simpler form of the Unear entropy Sr, it is 
easier to derive an approximate formula for its time evolution. Here we now present an 
analytical formaUsm for the time evolution of Sn in coupled strongly chaotic systems. Let 
us start the formahsm with the assumption that the initial state is a product state, given 
as IV'(O)) — ® |(/)2(0)), where |(/)i(0))'s are the states corresponding to individual 

subsystems. In general, the time evolution operator of a coupled system is of the form 
U = UeUo = Ue{Ui (g) U2), where Ue is the couphng time evolution operator and Ui's are 
the time evolution operators of the individual subsystems. Furthermore, we have assumed 
U, = cxp(— ieJfi2) where H12 = hM") ® h^'^\ and ft,^'^ are Hermitian local operators. 
Here we derive our formalism in the eigenbasis of /i*^*''s, i.e., /i^*) |ea^) = ei*'' |ea''), where 
{ea \ \ ea^)} are the eigenvalues and the corresponding eigenvectors of h^^K 

The one step operation of U on [V'(O)) will give the time evolving state at time n = 1, 
i.e., IV'(l))- Now, at this time, we can determine the matrix elements of RDM correspond- 
ing to the first subsystem in the eigenbasis of /i^'^ 's as 

[pi{l)U = E^^P (^«^ - '^'^] (ei'^ef |V'o(l))(V'o(l)|e^'\ef ), (22) 



where |V'o(l)) is the time evolving state of the uncoupled system. We now assume that 
IV'o(l)) is a random vector, and consequently we can further assume that the components 
of I ^/iQ ( 1 ) ) are uncorrected to the exponential term coming due to the coupling . Therefore, 
we have 



= (2) 



(23) 



where N is the Hilbert space dimension of the first subsystem and pio is the density matrix 
corresponding to the uncoupled top. If we follow same procedure for one more time step, 
then at n = 2 we have 



[Pi(2)] 



a/3 



(24a) 



where 
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Figure 5. Evolution of the Linear entropy for the coupled strongly chaotic system is 
presented. The dotted line is the numerical results of the coupled kicked tops system. 
We choose k = 6.0 for the first top and fc = 6.1 for the second top. The solid line is 
the theoretical estimation, given by Eq. 



a,/3 



,(2) 



(24b) 



If we evolve the initial state for any arbitrary time n and follow the similar procedure as 
described above, we get 



= (2) 



(25) 



It is now straightforward to calculate linear entropy from the above expression, and that is 
given as 

a,/3 7,(5 

For the coupled kicked tops H12 = Jzi ® Jz2 1 3- Therefore, for this particular system, the 
above formula would become in large j-limit as 



5K(n)~l-p(e)^("-i 



2 f, , Si(2iVe) ] / 1 
N\ e I V iVe 



X {1 - cos(2Are) + Ci(27Ve) - ln(27Ve) - 7} 



(27a) 
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where 



2 r 1 . [Ne 



) 



(27b) 



The functions Si and Ci are the standard sine-integral and cosine-integral function, respec- 
tively, while 7 = 0.577216... is the Euler constant. In the above formulation we have 
not assumed, unlike the perturbation theory [18], any particular order of magnitude of the 
coupling strength e. 

In Fig. 5, we have presented our numerical result of the linear entropy {Sr) production in 
the coupled tops where the individual tops are strongly chaotic. In the above formalism, we 
have not assumed any special symmetry property. Here we break permutation symmetry 
by taking two nonidentical tops with k = 6.0 for the first top and A; = 6.1 for the second. 
Fig. 5 demonstrates that how well our theoretical estimation, denoted by the sohd curve, is 
valid for both weak and strong coupling strengths. 



5. Summary 

In this paper, our major goal was to study entanglement production in coupled chaotic 
system as a possible quantum indicator of classical chaos. We have used coupled kicked 
top as a model for our study. Single kicked top is a well studied model of both classical 
and quantum chaotic system. Therefore, it is easier for us to identify the effect of under- 
lying classical dynamics on the entanglement production. We have studied entanglement 
production for different underlying classical dynamics of the individual top and also for 
different couphng strengths. Here we not only have considered entanglement production 
for the pure initial state, but we have also initiated the study of entanglement production 
in coupled chaotic systems for the mixed initial state. We have used Log-negativity, a re- 
cently proposed measure, to characterize mixed state entanglement production. In general, 
for both kind of initial states, entanglement production is higher for stronger chaotic cases. 
However, we have observed a saturation of entanglement production when individual tops 
are strongly chaotic and they are coupled strongly to each other. Coupling strength between 
two tops is also a crucial parameter for the entanglement production. For instance, when 
the coupling strength between two tops is very weak, we find higher entanglement produc- 
tion for sufficiently long time corresponding to non-chaotic cases. This is also a common 
property of both kind of initial states. We analytically estimated the above mentioned 
statistical bound on the pure state entanglement (saturation of entanglement production) 
using RMT We have also derived an approximate formula, based again on the ideas of 
RMT, for the pure state entanglement production in coupled strongly chaotic system. This 
formula is applicable, unhke perturbation theory, to large couphng strengths and is also 
valid for sufficiently long time. We still do not have deeper analytical understanding of 
mixed state entanglement production in chaotic systems, an open problem that is related to 
the establishment of calculable measures of mixed state entanglement. 
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